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ABSTRACT 



Context. Westerlund 1 is a young, massive Galactic starburst cluster that contains a rich coeval population of Wolf- 
Rayet stars, hot- and cool-phase transitional supergiants, and a magnetar. 

Aims. We use spectroscopic and photometric observations of the eclipsing double-lined binary WIS to derive dy- 
namical masses for the two components, in order to determine limits for the progenitor masses of the magnetar 
CXOU J164710. 2-455216 and the population of evolved stars in Wdl. 

Methods. We use eleven epochs of high-resolution VLT/FLAMES spectroscopy to construct a radial velocity curve for 
WIS. 7?-band photometry is used to constrain the inclination of the system. 

Results. WIS has an orbital period of 9.2709±0.0015 days and near-contact configuration. The shallow photometric 
eclipse rules out an inclination greater than 65°, leading to lower limits for the masses of the emission-line optical 
primary and supergiant optical secondary of 21.4 ± 2.6Mq and 32.8 ± 4.0Mq respectively, rising to 2S.2j;3 QMQ and 
35.41|;°Mq for our best-fit inclination 62I4 degrees. Comparison with theoretical models of Wolf-Rayet binary evo- 
lution suggest the emission-line object had an initial mass in excess of ~S5M0, with the most likely model featuring 
highly non-conservative late-Case A/Case B mass transfer and an initial mass in excess of 40Mq. 

Conclusions. This result confirms the high progenitor mass of the magnetar CXOU J164710. 2-455216 inferred from its 
membership in Wdl, and represents the first dynamical constraint on the progenitor mass of any magnetar. The red 
supergiants in Wdl must have similar progenitor masses to WIS and are therefore amongst the most massive stars to 
undergo a red supergiant phase, representing a challenge for population models that suggest stars in this mass range 
end their redwards evolution as yellow hypergiants. 

Key words, stars: evolution - supergiants - stars; individual: WIS - stars: magnetars - binaries: general 



1. Introduction 

The Galactic starburst cluster Westerlun d 1 (hereafter 
Wdl: I Westerlund lll96lUClark et al~ll2005D contains a rich 
coeval population of OB supergiants, yehow hypergiants 
(YHGs) and red supergiants (RSGs) that collectively map 
out the transitional post-Main Sequence (MS) loop red- 
wards followed by massive stars in the cluster. The relative 
brevity of the transitional phase and intrinsic rarity of such 
objects means that this stage of evolution is poorly under- 
stood, but recent downward revisions to MS mass loss rates 
(|Fullerton at al. II2006I: [Mokiem et al. II2007D suggest that it 
is of critical importance in understanding how massive stars 
shed their outer layers prior to the Wolf-Rayet (WR) phase. 

Due to its unique stellar population, Wdl has been 
the subject of intensive observational study in recent 
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years (jRitchie et ah 2009a, hereafter Paper I, and refs. 
there in: see also IClark et al. I l2010al iNeeueruela et al. 1 
I20T0I) . Studies of the massive stellar population support 
a single burst of star formation at an age ^5Myr and a 
distance ~5kpc ijCrowther et al. 1 120061 : INeeueruela et al. 1 
i20100 . with the identification of '^08V stars in the clus- 
ter (Clark et al., in prep.) and a population of lower- 
luminosity late-0 II-III stars just evolving off the MS both 
fully consistent with this derive d age. Dynamical mass de- 
terminations of late-0 dwarfs (|Gies I [2003i) and compar- 
ison of the popul ation of OB stars i n Wd l with theo- 
retical isochrones (jMevnet fc Maeder 1 12000|) suggest that 
stars with Mini~3OM0 lie at the MS turn-off, with the 
early-B supergiants having pro genitor masses Mini~35- 
4OM0 (|Negueruela et al. I I2OIOI ) and the WR population 
desce nded from stars with Mini>4OM0 (jCrowther et al. I 
120061 ). However, to date no direct mass determination exists 
for a member of Wdl; as well as providing confirmation of 
the current understanding of the cluster derived from spec- 
troscopic studies, this is of importance for confirming the 
high progenitor mass for the magnetar CXOU J164710.2- 
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Fig. 1. Left panel: /-band spectra of W13 at (p — 0.33 and 
4> — 0.74. Right Panel: phase dependent variability of the 
Pal2 emission line. Five spectra are overplotted, with three 
from (j) = 0.15 — 0.30 (solid line) and two from (j> — 0.30 — 
0.35 (dotted line). 



455216 that has be en inferred from its membership of Wdl 
uno et ai71l2006[ ). In addition, the distribution of evolved 
stars in Wdl offers the prospect of demanding tests for evo- 
lutionary models, with both the distribution of WR sub- 
types and the large number of cool hypcrgiants in both 
YHG and RSG ph ases at odds with current predictions 
(|Clark et al. l[2010al) . 

In this paper we present spectroscopic radial velocity 
(RV) measurements of the mas sive binary W13 , identified 
as a 9.2-day eclipsing system bv lBonanos I (|2007[ ). X-ray ob- 
servations reveal a hard source wit h L^j^lO'^^ergs"^ consis- 
tent with a colliding- wind system (jClark et al. 1120081 ) . while 
subsequent multi-epoch spectroscopy described in Paper I 
showed W13 to be a double-lined spectroscopic binary con- 
sisting of a B0.5Ia+/WNVL emission-line objeclQ and an 
early-B supergiant. We supplement the results of Paper I 
with an additional six epochs of data, providing a total 
baseline of approximately 14 months that allows an accu- 
rate radial velocity curve to be derived for the W13 sys- 
tem. This is use d in conjunction with i?-band photometry 
(jBonanos Il2007f ) to determine the masses of the two com- 
ponents. 



ratio of our spectra is ~95 at 8700A . Photome t ry of W13 
was taken from the published data of 'Bona nos I (|2007(l , ob- 
tained on 17 nights between 15/6/2006 and 25/7/2006 us- 
ing the Im Swope telescope at Las Campanas Observatory, 
Chile. 

RV measurements were carried out using the IRAlfl 
onedspec tasks to fit Gaussian profiles to the cores of the 
Paschen series absorption and emission features in the spec- 
trum of W13: the derived RV at each epoch is an error- 
weighted average of the RV of individual lines. Strong in- 
terstellar C2 lines from the (2,0) Phillips band overlap the 
Pal2 line, leading to a systematic phase-dependent bias of 
up to ^lOkms^"'^ in the line centre, and this line was there- 
fore excluded from t he analysis. In addition, a broad DIB 
centred at ~ 8648 A ()Negueruela et al.1l20lOD overlaps the 
Pal3 A8665 line, attenuating the blue flank of the emis- 
sion component at (^'^0.25 and leading to an offset red- 
wards relative to the other Paschen series lines, although 
measurement of the absorption component is not affected. 
Finally, at t/)'^ 0.1 5-0. 3 a systematic decrease in strength is 
seen in all Paschen-series emission lines, with the lines re- 
covering in strength rapidly after 0~O.3; this can be seen 
in the right panel of Figure [U which overplots five spec- 
tra taken at 0~0. 15-0.35. The weakening in Paschen-series 
emission was assumed to be wind variability in Paper I, but 
the persistence of this behaviour over more than ~30 orbits 
in our extended dataset implies that it is not a transitory 
effect but rather a region of excess absorption periodically 
crossing the line of sight. 

Absorption line RV measurements were therefore car- 
ried using the Pall and Pal3-16 lines, noting that the 
progressive weakening of the higher Paschen-series lines 
leads to an increased fitting error and consequent decreased 
weight in the derived RV. Emission line RV measurements 
were carried out as follows: 

— Emission lines cannot be measured with accuracy for 
one spectrum taken near eclipse (0=0.03, 24/07/2008, 
MJD=54671.13). 

— At 0=0.1-0.3 attenuation of the blue flank of Pal3 by 
the ~8648A DIB and the decrease in strength of Pal 5 
and Pal6 leave only Pall and Pal4 emission lines avail- 
able for measurement, with the two lines in good agree- 
ment. 

— At other epochs. Pall and Pal3-16 emission compo- 
nents were measured. 



2. Observations & data reduction 

Spectra of W13 were obtained on eight epochs in 2008 and 
three epochs in 2009 using the Fi bre Large Array Mul ti 
Element Spectrograph CFLAMES: IPasauini et aTl 1200211 . 
located on VLT UT2 Kueyen at Cerro Paranal, Chile. The 
GIRAFFE spectrograph was used in MEDUSA mode with 
setup IIR21 to cover the 8484-9001A range with resolution 
i?~ 16200; full details of data acquisition and reduction are 
given in Paper I, and representative spectra showing the 
Pall A8862 and Pal2 A8750 lines at two extremes of the 
RV curve are shown in Figure [TJ The signal-to- noise {S / N) 

^ The WNVL classification is used here to indicate strong 
spectroscopic similarities with other WNL stars in Wdl; we 
note the He II A4686 lin e necessary for formal classification 
l|Crowther fc Smith 11199^ ) lies outside our spectral coverage. 



3. Results 

3.1. Spectroscopic classification 
3.1.1. The optical primary 

The i?-band spectrum of W13 is dominated by 
the emission-line optical primary, with strong, rela- 
tively broad Ha emission (FWHM'^500kms~'^), com- 
plex, time-varying P Cygni profiles extending to at 
least — 350kms~^ in the He I AA6678, 7065 lines and 



^ IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooperative 
agreement with the National Science Foundation. 
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Table 1. Journal of observations. 



Date 


MJD" 


Phase* 


Elapsed orbits* 


RV„h„ fkms"^l 


RVcm (kms"M 


on ICiC^ /onnQ 
zy/UD/zUUo 


KARA a iQ/if; 
04D4D.io4D 


U.oo 


u.oo 


1 rr /j 1 n 

+o4±y 


non 1 ^ 
— ZOZ± / 


1 o /AT /onno 
io/U < /2UUe 


o4ddo.Uood 


U.O ( 


2.3 ( 


1 on 1 o 


— 24D±io 


24/07/2008 


54671.1343 


0.03 


3.03 


-50±11 




14/08/2008 


54692.0423 


0.28 


5.28 


+74±13 


-224±14 


04/09/2008 


54713.0107 


0.55 


7.55 


-120±5 


+ 127±29'= 


15/09/2008 


54724.0818 


0.74 


8.74 


-200±9 


+163±11 


19/09/2008 


54728.0554 


0.17 


9.17 


+56±4 


-225±17 


25/09/2008 


54734.0613 


0.82 


9.82 


-173±7 


+164±9 


14/05/2009 


54965.1768 


0.76 


34.76 


-214±14 


+137±8 


18/05/2009 


54969.3198 


0.21 


35.21 


+76±13 


-242±7 


20/08/2009 


55063.0575 


0.32 


45.32 


+78±11 


-217±16 



'"Modified Julian day at the midpoint of two 600s integrations (2x500s, 04/09/2008; Ix600s+lx700s, 19/09/2008). 
'Phase taking ro=54643.080, elapsed orbits for an orbital period of 9.271d (see Section [3. 2|) . 

'^Unreliable due to the appearance of a central reversal in the emission line, and excluded from the fit to the RV curve. 
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Fig. 2. Intermediate-resolution spectra of W13 and the blue 
hypergiants W5 (B0.5Ia+/WN1 0-ll, top) an d W7 (B 5Ia+, 
bottom). Spectra taken from IClark et al. I (|2010al ) and 
iNeeueruela et al. I (|2010f ). 



weak C II AA6578,6582 emission. A comparison of pre- 
viously published intermediate-resolution i?-band spec- 
tra of W13 and the blue hypergian ts W5/WR S 
('B0.5Ia +/WN10-ll: INeeueruela fc Clark I l200g) and W7 
fBSIa"*": [Neeueruela et al 2010ll is shown in Figure [21 The 
luminous mid-late B hypergiants in Wdl EI are all believe d 
to be on a pre-RSG loop redwards ( Cla rk et al. Il2010b() . 
and display narrow Ha profiles with strong P Cygni ab- 
sorpt ion components, along with He I and C II in absorp- 
tion (jClark et al. Il2010al: INeeueruela et al. Il2010t ). In con- 
trast, the spectrum of W13 shows strong similarities to both 
W5 and the WN9h bina ry W44/WR L ( Crowthc r et al. I 
l2006tlcrark et al■1l20l0a^ ■ with the three objects forming a 
morphologically-distinct group when compared to the mid- 
late B hypergiants, suggesting that the emission-line object 
in W13 is also an immediate evolutionary precursor to the 
WR phase. 



The spectrum plotted in Figure [5] shows He I lines that 
are redshifted by ~50kms~^ and C II emission lines some- 
what bluewards of He I (RV^Okms^^), although accurate 
measurement of the C II line centres is difficult due to 
their low strength and the relatively low S/N of the R- 
band spectrum. Nevertheless, if the C II lines originated in 
the companion then we would expect the lines to be sig- 
nificantly blueshifted at this epoch (RV^— 150kms~^; see 
Section [3.2p . suggesting that the emission lines have a com- 
mon origirO- However, N II A6611 emission, present in both 
W5 and W44, appears absent in W13, while the Ha, He I 
and C II emission lines are also considerably weaker, sug- 
gesting W13 is the least evolved member of the WNVL 
population of Wdl. 



3.1.2. The optical secondary 

The Paschen-series lines display complex emis- 
sion/absorption profiles with the two components moving 
in anti-phase (see Figure [T]). The absorption components 
have similar strengths to the lower-luminosity 09.5Iab/b 
objects in our FLAMES dataset, but infilling from the 
emission-line object is likely to affect these features and 
an alternative diagnostic is provided by weak He I A8584, 
8777 absorption lines that are apparent in many spectra, 
moving in phase with the Paschen-series absorption lines. 
These He I lines ar e first seen at 08-91 an d strengthen 
rapidly at -61.51 (jNeeueruela et al. I l20100 . with their 
weakness therefore suggesting a spectral type no later than 
~B1I, although overlapping interstellar features (including 
a broad, weak DIB at ~8779A) preclude precise mea- 
surement. A late-0 spectral type appears to be excluded 
by the apparent absence of C HI A8500 absorption, with 
this line leading to a bluewards offset and discrepant 
strength for the C HI A8500/P al6 A8502 blend i n the 
O9.5-B0.5 supergiants in Wdl (Neeueruela et al. I I20T0I: 
see also Paper I). Neither classifier should be significantly 
affected by wind emissior^, and the weak He I lines and 



^ We note that W5 also appears to show C II blueshifted rel- 
ative to He I. 

^ The He I A8584,8777 lines arise from transitions from the 

W7, W33, and W42a; B5-B9Ia+ (|Negueruela et al. l[20Tol '). 3p^P° level, which is well-populated in ~Bl-4 supergiants, to 
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Fig. 3. Radial velocity curve for W13. 



Fig. 4. Fit to the R-band lightcurve of W13 



absence of C III absorption therefore suggests a ~B0.5-1I 
classification with an uncertainty of roughly half a spectral 
subtype and lack of strong constraints on the luminosity 
class. However, we caution that if C III is weak due 
to abundance anomalies or near-critical rotation then 
this limit may not apply, and the weakness of the He I 
lines permit a classification of 09.5-BOI that is broadly 
consistent with the Paschen-series line strengths. 

3.2. Radial velocity curve 

Figure [3] shows RV curves for the two compo- 
nents of the W 13 systein . Taki ng the 9.20-day pe- 
riod reported by iBonanos! (|2007t) as a starting point, 
an error-weighted fit to the radial velocities of 
the absorption-line component yielded best-fit val- 
ues for the orbital period of P = 9.2709 ± 0.0015 
days, consistent with an independe nt determination 
using a Lomb-Scargle periodogram (j Press fc Rvbicki I 
Il989|) . a systemic velocity of —65.9 ± 2.4 kms~^ and 
semi-amplitude i^abs = 137.3 ± 6.7 kms~^. The corre- 
sponding fit to the emission line RV curveQ has 
a systemic velocity of —48.2 ± 3.1 kms^^ and semi- 
amplitude it'om = 210.2 ± 8.7 kms^^. Errors are derived 
from the fitting residuals using the bootstrap method 
(|Efrom fc Tibshirani 11199"^ . We note that systemic veloc- 
ity derived from the emission line fit is somewhat lower 
than that derived from the companion, and is in closer 
agreement with the mean radial velocities of other Wdl 
supergiants (see Paper I). Discrepancies in this parameter 
are commonly observed in early-type spectroscopic bina- 
ries (e.g. the 9.8- day 07III(f)-H08.5I binary HD 149404; 
iRauw et al. Il200l] ) although the effect is small in compar- 
ison with some other e volved systems (e. g. HDE 228766; 
iMassev fc Contil 119771 : iRauw et al."1 12002[ ) in which wind 
contamination strongly affects derived systemic velocities. 

Taking these values yields a mass ratio 
g = Mabs/A/cm = 1-53 ± 0.10 and masses for the two 



components of: 



upper levels near the ionisation limit of He I ()van Helden \\l97^ ). 
but the corresponding downwards transitions are very weak. 

® In fitting the emission-hne RV curve, the discrepant point 
at <h=0.55 was excluded. 



sin^i 



27rG 



15.9± 1.9Mt 







and 



Mabs sin'^z = 24.4 ± 3.0Me 



(1) 



(2) 



Finally, measurements of blended hydrogen lines with 
Gaussian fits tend to yield systematically lower values of 
Ki and K2 than meth ods such as spectral disentangling 
(ISimon fc Sturm _ 1994) that are less affected by blend- 



ing ()Andersen I Il975t ISouthworth fc Clausen I [20071 ). The 

paucity of strong lines free from significant interstellar, tel- 
luric and wind contamination in the R- and /-band spec- 
tra of ~B0 supergiants makes the extent of this effect on 
our determination of Kcm and i^Tabs hard to quantify, and 
we therefore note that our method may underestimate the 
masses of the two components of W13. 

3.3. Light curve 

To constrai n sin^z, we folded the i?-band photometric data 
reported bv lBonanos I (|2007t ) on to the 9.271 day period de- 
termined from the RV data. The data were binned to reduce 
the considerable scatter present in the light curve, which is 
most probably a consequence of intrinsic aperiodic vari- 
ability in one or both components: low-level photometric 
and spectroscopic variability is a feature of all transitional 
supergiants in Wdl, with the blue hypergiants displaying 
rapid photometric variability at the '^O.l magnitude level 
and the early-B sup ergiants also variable (jBonanos 1 120071 : 
IClark et al. Il2010al ). Therefore, given the limited dataset 
and shallow ~0.15 magnitude eclipse, we do not expect to 
be able to fit the light curve of W13 with high accuracy. 
Nevertheless, in the ab sence of l o nger- term photometric 
monitoring, the data of iBonanos I (|2007[ ) allow reasonable 
constraints to be placed on the orbital inclination. 

The nightfall cod^II was used to model the light curve 
of W13. The effective temperature of the emission line ob- 
ject was fixed at 25kK, appropriate for its spectral type, 
and a linear limb-darkening law and circular orbits were as- 
sumed. The mass ratio was derived from the RV curve, and 



http:/ /www.hs.uni-hamburg.dc/DE/Ins/Pcr/Wichmann/Nightfall 
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the indination, Roche lobe filling factors for both objects 
and temperature of the optical secondary were allowed to 
vary. The code rapidly converges to a near-contact configu- 
ration in which the emission-line object has almost filled its 
Roche lobe (filling factor ~0.93±0.05) and the other star 
has a somewhat lower filling factor (~0.74±0.1). The best- 
fit model has an inclination i — 62° and provides a close 
match to the light curve from (/)~0.5 to 0^1, although the 
region around (p = 0.2-0.4 is less well reproduced. Although 
the nightfall code supports additional features such as 'hot 
spots' that may provide a better fit to this portion of the 
light curve ( see, for example, th e model of Cyg OB2#5 
presented bv iLinder et al. I [20091 ) . we consider further re- 
finement of the model inappropriate given the limitations 
of the photometric dataset used. Parameters derived from 
the light curve model are listed in Table [H although we 
stress that our primary goal is to constrain the inclination 
of the system and other parameters should be regarded as 
provisional pending acquisition of longer-term photometry. 

To examine errors in the derived inclination we investi- 
gated models in which i is fixed at values from 55° to 68° 
while filling factors and temperatures are allowed to vary as 
before; the best fit i = 62° curve and models with i = 59°, 
i — 65° and i ~ 68° are plotted in Figure S) Inclinations 
greater than our best-fit model lead to slightly lower filling 
factors and higher temperatures, while the converse is true 
for lower inclinations. The depth of the eclipse provides 
the strongest constraint on the model, with inclinations 
greater than ~65° stro ngly disfavoured , /-band photom- 
etry is also presented bv lBonanos I ()2007t ). but considerably 
greater scatter is present in the data at mid-eclipse, making 
it less suitable for modelling. Nevertheless, these data also 
support i < 65°, favouring a value '^60-62°, although the 
degree of scatter renders this uncertain. Therefore, taking 
i < 65° places robust lower limits of Afcm > 21.4 ± 2.6M0 
and Mabs > 32.8 ± 4. CM© for the emission-line object and 
its companion, rising to 23.2j]3'QM0 and 35.4I4 gM© for our 
preferred inclination i — 62+^° . For the purposes of discus- 
sion we take the lower masses derived from the i < 65° 
limit. 

4. Discussion and Conclusions 

4.1. The evolution of the W13 system 

The short orbital period, near-contact configuration and 
evolved, mass-depleted nature of the emission-line object 
all imply that the two components of the W13 system 
must have undergone strong interaction during their evo- 
lution. The 9.3-day orbital period suggests a late-Case A 
or Case B scenario, with mass tra nsfer beginning near 
the onset of shell hydrogen burning (jPetrovic et al. 1120051 
hereafter P05). The presence of unevolved late-0 stars 
(Mini ~ SOMq) in Wdl suggests a minimum initial mass 
~35M0 for shell burning to have commenced, with in ex- 
cess of ~1OM0 lost once mass transfer begins. However, 
transfer of angular momentum is expected to lead to the 
accretor rapidly re aching critical rotation (jPacket I Il98ll : 
iLanger et al. II2008I) . while rap id rotation wi ll also greatly 
increase wind mass loss rates ()Langer II1998D . and a fully- 
conservative transfer scenario appears unlikely. Indeed, 
models of short-period WR-l-0 binaries by P05 suggests 
that mass-transfer is highly non-conservative in such sce- 
narios, with only ^-^10% of transferred mass being retained 



Table 2. Summary of orbital and physical parameters of 
W13''. 



Parameter 


Value 




iQ [NUD) 
(.daysj 


^ Ad AO non 

OTTio -1- n nni k 
y.z (Uv ± U.UUlo 




q ~ JTlabs /niom 
n ( R^) 

g 
i 


1 CO _|_ n 1 n 

/ zzco 
('fixed ^ 
62+1° 






Emission 


Absorption 


Filling factor 


0.93 ±0.05 


0.74 ±0.1 


Teff (K) 


25000 (fixed) 


25000±2000 


R (Re) 


22±2 


21±2 


7 (km s 


-48.2 ± 3.1 


-65.9 ±2.4 


K (km s"^) 


210.2 ±8.7 


137.3 ±6.7 


Msin^i (Mq) 


15.9 ±1.9 


24.4 ± 3.0 


M{i = 65°) (Mq) 
M{i = 62lf ) (Mq) 


21.4 ±2.6 

OQ 9+3.3 
^■^•^-3.0 


32.8 ±4.0 

35.411° 



" Note that To corresponds to the eclipse of the BO. 51a''' 
emission-line star. 



by companion star. The current 21Mq+33Mq (minimum) 
mass ratio in W13 is consistent with the P05 model of late- 
Case A/Case B evolution at low accretion efficiency, with 
higher accretion efficiencies leading to a more unequal mass 
ratio than we observe. 

P05 estimate a relationship between initial MS 
mass and final WR mass for Case B systems of 
Mini = (MwR + 4.2)/0.53 that suggests that the emission- 
line object had an initial mass ~48Mq. While this 
i s consistent with est imates of WR progenitor masses 
(|Crowther et al. Il2006[) . it is somewhat higher than ex- 
pected for a star just entering the WR phase in Wdl. 
However, W13 will likely shed further mass before becom- 
ing a bona fide WN9, suggesting that Mwr ^ 20Mq and 
consequent initial MS mass around ~45M0 are more ap- 
propriate. Assuming an accretion efficiency of ~10% from 
P05, this would imply that ~2-3M0 was transferred to 
the secondary, with the remainder lost from the system. 
While extended radio e mission from W13 is not detected 
(jPoughertv et al. Il2010 t). emission from both the 09Ib star 
W15 and the extreme RSG W26 overlaps the region around 
W13 and might obscure direct signs of recent mass loss. 

W13 therefore appears to be a less-evolved analogue to 
WR21 (HD 90657), a 8.3-day 19±37Mo WN5 binary con- 
sidered by P05. One notable discrepancy is the presence 
of an evolved companion in W13, whereas WR21 contains 
an unevolved mid-0 star. Although abundance anomalies 
might suppress the C HI A8500 line in a late-0 supergiant, 
leading to an erroneously late spectral type, the presence 
of He I lines moving in phase with the Paschen-series ab- 
sorption lines excludes a spectral type earlier than ~09 
(jNeeueruela et al"~ll2010l) . It is possible that the supergiant 
is still in extreme, near-critical rotation from recent mass 
transfer, and is thus expanded with a complex, latitude- 
dependent spectrum that features both hot (polar) and 
cooler (equatorial) components. However, confirmation of 
this hypothesis is observationally challenging, as infilling of 
the Paschen-series absorption lines prevents direct deter- 
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Fig. 5. Location of luminous hypergiants in Wdl for a 
distance of 5kpc compared to z=0.020 (i.e. solar abun- 
dance) evolutionary tracks with and without rotation 
(|Mevnet &: Maeder I [2003). Lower limits for the luminosi- 
ties of the YHGs W4, W16a and W265 are plotted. W4 
is offset vertically by 0.02 dex and error bars are omitted 
for W265 to highlight the change in spectral type over the 
pulsational cycle. 



mination of v sini while the high reddening towards Wdl 
precludes the use of fiducial 0-type spectral classifiers in 
the blue region of the spectrum. 

4.2. Evolutionary implications for Wdl 

4.2.1. The cool hypergiant population of Wdl 

The likely >4OM0 MS mass for the emission-line object 
in W13 provides direct constraints on the masses of the 
eleven cool hypergiants in Wdl. Of these, only two have 
be en studied in d etaiH but the long-term dataset compiled 
by IClark et al. I (|2010al) reveal a remarkable lack of secu- 
lar evo lution amongst these objects, with only the LBV 
W243 (|Ritchie et al. Il2009b[) apparently under going a ma- 
jor outbu rst in the half-centu ry since the discovery of the 
cluster bv I Wester lund I (|1961[ ). Although early observations 
are sparse, the available data are sensitive to the long- 
term evolutiona ry trends seen in objects such as the YHG 
IRC -HO 420 (iHumphrevs et al. I [200l. the LBV R,127 
(IWalborn et al. Il2008[ ) or M33's Var A (jHumphrevs et al. I 
I2006D . Current observations therefore suggest the super- 
giants in Wdl undergo a slow redwards evolution at ap- 
proximately constant luminosity until they encounter an 
extended cool-ph ase state accompani ed by growing pulsa- 
tional instability (jClark et al. Il2010allbh . 

Evolutionary models of massive stars predict a rela- 
tively long YHG lifetime at T ~ BOOOX for stars with 
Mini~40M(7), consistent with the lack of secular evolution 



and mid- A to late-F spectral types for the YHGs in Wdl, 
but do not predict further evolution to the RSG phase 
(jMevnet fc Maeder I 120031: iDrout et al.1 I2009D . However, 
while the B- and A/F- hypergiant populations plotted in 



FigureJS] appear in good agreement with the evolutionary 
track^, the presence of RSGs in Wdl suggests that stars 
in this mass range do evolve further redwards (Clark et a i~| 
l2010a), with the extended ra dio nebulae around these ob- 
jects ( Doughertv et al. ||2010D revealing extensive mass loss 
that may account for a significant fraction of the mass 
lost prior to the WR phase. Lower-mass stars will not 
reach the RSG phase until >6Myr, and the absence of 
any other indicators of significant non-coevality in Wdl 
(jNegueruela et al.1 120101 ) and the location of two RSGs 
near the core of the cluster argue against the RSGs be- 
ing descended from a separate population of older, lower- 
mass stars. Short-term spectroscopic variability and lack 
of contemporaneous photometry render the luminosities 
of individual RSGs uncertain, and further observations 
are required to place them firmly on the HR diagram. 
Nevertheless, the Ml-5Ia spe ctral types derived f rom TiO 
band strengths appear robust (jClark et al. Il2010"al) . and are 
clearly discrepant with respect to current theoretical pre- 
dictions. 



4.2.2. The Wdl magnetar. 



The magn etar CXOU J164710.2-455216 (|Muno et al. I 
I2006L I2007D lies 1'.7 from the centre of Wdl, correspond- 
ing to 2.3((i/5kpc) pc in projection. The negligible likeli- 
hood of the magnetar being a chance association and the 
presence of ~08V stars in Wdl prov ide strong evidence 
for a massive progenitor uno et al. 112006 ). a resuh con- 
firmed by our direct measurement of the mass of W13 
which rules out a magnetar progenitor below '^35M0 and 
strongly supports a progenitor mass in excess of ~45Mq 
unless mass transfer within W13 is unexpectedly conser- 
vative. This is consistent with the expected initial masses 
of the most evolved WR stars in Wdl, and also pro- 
genitor masses derived for other magnetars: a 48^^ M© 
progenitor mass is found for the magnetar SGR1806-20 
based on its assumed members hip of a massive cluster 
at GlO.0-0.3 (jBibbv et al. II2008D . while an expanding HI 
shell around the magnetar IE 1048.1-5937 is inferred to 
be a wind-blown bubble from a ~3O-4OM0 progenitor 
(Gaensler et al. 2005). Th e progenitor o f SGR19 00+14 has 
a sign ificantly lower mass (IClark et al. ll2008l: lD avies et al. I 
l2009f ). implying a number of formation pathways exist. 
However, the derived mass of SGR1806-20 is strongly in- 
fiuenced by the distance to the host cluster obtained from 
spectroscopic observations, and our results place the first 
dynamical constraints on a massive magnetar progenitor. 
Indeed, our resuhs suggest that CXOU J164710. 2-455216 
may have the highest dynamically-constrained progenitor 
mass of any confirmed neutron star. The high-mass X- 
ray binary 4U1700-37 contains a 06.5Iaf+ mass donor with 
mass 58±11Mq but the nature of the 2.44±q.27M com- 
pact object remains uncertain (jClark et al. I l2002l ). while 
the ~43±1OM0, highly-luminous Blla+ hypergiant in the 
X-ray binary Wray 977/GX 302-2 may have evolved via 
quasi-conservative Case A transfer in a 3.5 day binary with 



* The A3Ia+ LBV W243 ([Ritchie et al. ll2009bD and the YHG 
W265, which varies from Fl-5Ia+ with a -100 day quasi-period 
(Paper I; also IClark et al. Il2010al ') 



^ Luminosities f or the BHGs W7 and W33 are taken from 
iNegueruela et al. I ([2010), while the luminosity of the LBV W243 
is taken from non-LT E modelling assuming a distance of 5kpc 
(jRitchie et al. Il2009bl '). YHG luminosities of log(L/L0) = 5.7 
are lower limits from the My-W(0 I A77 74) relationship 
(jArellano Ferro et al. |[2003 : IClark et al.1l2005D . 
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initial masses 26M0+25M0 (IWellstein fc Langer I 119991: 
iKaper et al~ll2006l^ . 

NTT/SofI i^'s-band ima ging rules out a cur rent > IMq 
companion to the magnetar (jMuno et al. I[2"006l ). a result we 
confirm using deep VLT/NACO J-, H- and i^s-band imag- 
ing (Clark et al., in prep.). Nevertheless, g iven the high bi- 
nary fraction amongst the WR population ( Crow ther et al. I 
[200a) and the need for a low pre-supernova core mass 
to avoid direct ( or fallback) black hole formation (e.g. 
iFrver et al. 1 120021 ) it would appear likely that the mag- 
netar progen itor was part of a (now-disrupted) close bi- 
nary system (jClark et al. Il2008f ). Support for this hypothe- 
sis comes from population synthesis models, which can only 
form a neutron star from an isolated ~6OM0 progenitor 
within the ~5Myr age of W dl if mass loss ra tes from stellar 
winds are greatly enhanced (jBelczvnski fc T aam 2008). In 
a close binary scenario, however, removal of the hydrogen- 
rich outer mantle via Case A mass transfer results in a 
reduced post-MS helium core, and ongoing Case B trans- 
fer during shell burning will leave a low mass (< IOMq) 
helium-burning WR (P05), permitting isolated neutron star 
formation within 5Myr via a type Ib/c supernova if the kick 
velocity is sufficient to disrupt the svsterrFl. Neutron star 
formation may also occur for massive binaries with initial 
periods greater than a few weeks; such systems will not 
undergo Roche lobe overflow until core hydrogen burning 
is complete and therefore form higher-mass helium cores 
than Case A systems, but if Case B or early-Case C mass 
transfer can suppress hydrogen shell burning before core 
helium burning is complete then the consequent reduction 
in th e mass of the iron core may limit black hole forma- 
tion (|Brown et al. |[200l[ ). Further results on the distribu- 
tion of binaries on Wdl from our VLT/FLAMES survey 
and follow-up observations will therefore allow strong con- 
straints to be placed on the formation channels for such 
systems. 

5. Conclusions and future work 

We find lower mass limits for the components of the eclips- 
ing binary W13 of 21.4 ± 2.6Mo and 32.8 ± 4.OM0, ris- 
ing to 23.2t^;j^MQ and 35.4t^;°MQ for our best-fit incli- 
nation 62^4 degrees, with spectroscopy suggesting that 
the evolved emission-line object is likely an immediate 
evolutionary precursor to the WR phase. As conservative 
mass transfer would require the exchange of (at least) 5- 
IOMq without the accretor exceeding critical rotation, it 
appears likely that W13 evolved through non-conservative 
late-Case A/Case B mass transfer as the (initially) more 
massive star left the MS. Estimates of the initial WR 
mass from P05 and the presence of a >33M0 super- 
giant companion, which cannot have greatly increased in 
mass during highly non-conservative transfer, therefore sug- 
gest a MS mass for the emission-line object in excess 
of ~4OM0. This is consistent with previous estimates of 
the transitional supergiant masses in Wdl obtained from 
the MS turnoff and s pectroscopy of the WR and OB su- 
pergiant populations (ICrowther et al. I l2006t IClark et al. I 
l2010at iNegueruela et al. 1 12010'). Most importantly, this re- 
sult places the first dynamical constraint on the mass of a 

In this scenario, ~90% of the MS mass of the primary is lost 
prior to SN, with the second ary ultimately forming an isolated 
black hole with Mbh~8Mq (|Belczvnski fc Taam |[200g 'l. 



magnetar progenitor, and highlights a discrepancy between 
the presence of RSGs in Wdl and the predictions of evo- 
lutionary models, which suggest that the most luminous 
RSGs should evolve from significantly lower masses. 

A first study of the binary fraction amongst lower- 
luminosity late-0 II-III stars in Wdl will be presented in 
a subsequent paper in this series, but many binary sys- 
tems are already available for follow-up study. These in- 
clude short-period spectroscopic binaries (W43a, W3003; 
Paper I) , eclipsing binaries withi n the WR, OB s upergiant 
and main sequence populations ()Bonanos 1 12007[). and X- 
ray and radio-selected coUiding-wind binaries (jciark et al. I 
.2008; .Dougherty et al. ,,2010.) . Consideration of these data 
will allow further dynamical constraints to be placed on the 
progenitor masses of the evolved stars within Wdl as well 
as the general mass luminosity relation for stars in the up- 
per reaches of the HR diagram and the post-MS pathways 
they follow. Moreover, they will yield the first characterisa- 
tion of the binary properties of a homogeneous population 
of massive stars, of critical importance for studies of both 
star and cluster formation and numerous high-energy phe- 
nomena such as supernovae. Gamma-ray bursters and the 
formation of high mass X-ray binaries. 
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